Dynamics of Anderson localization in open 3D media 
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We develop a self-consistent theoretical approach to the dynamics of Anderson localization in 
open three-dimensional (3D) disordered media. The approach allows us to study time-dependent 
transmission and reflection, and the distribution of decay rates of quasi-modes of 3D disordered 
slabs near the Anderson mobility edge. 



Because of scattering from heterogeneities, the trans- 
mission of a short wave pulse through a disordered 
medium is strongly delayed. The precise shape of 
the decay of the average transmission coefficient T(t) 
has been extensively studied during the recent years 
IH 0. El Eli parallel to related studies of current 
relaxation in disordered conductors 0, El Ell E3 and 
survival probability decay in classically chaotic open dy- 
namical systems |l2j . If the disorder is weak, wave prop- 
agation is diffusive and T(t) decays exponentially with 
time t 0,0- This simple result only holds, however, un- 
til the Heisenberg time tji — 1/A, where A is the typical 
inter-mode spacing. Beyond tjj the decay of T(t) is be- 
lieved to be slower than exponential due to the so-called 
"pre-localized" states (or modes) that have anomalously 
long lifetimes [jj, 13, 111 |l(J, I 111. For strong disorder, An- 
derson localization sets in [lA. H3. flil] and the exponen- 
tial decay is expected to disappear, at least for quasi- 
one-dimensional (quasi-lD) samples where tu becomes 
smaller than the arrival time of the pulse. Much less is 
known about the effect of Anderson localization on T(t) 
in three-dimensional (3D) samples, where the Ioffe-Regel 
criterion of localization k£ < 1 is satisfied or, at least, 
approached as in recent experiments by Johnson et al. 
|6j (k is the wavenumber, £ is the mean free path due 
to disorder). None of the existing analytic approaches to 
Anderson localization (random matrix theory [l6j |. non- 
linear cr-model 0,E10]> optimal fluctuation method [Icj . 
self-consistent diagrammatic theory have succeeded 
so far to model the dynamics of Anderson localization in 
open 3D media. 

In the present Letter we apply an improved version of 
the original self-consistent theory of localization by Voll- 
hardt and Wolfle to calculate the transmission and 
reflection coefficients, as well as the distribution of mode 
decay rates of open 3D disordered media in the localized 
regime. The basic idea in this theory is to renormalize 
the diffusion coefficient to account for constructive inter- 
ferences between reciprocal wave paths inside the sam- 
ple. The self-consistent theory of Anderson localization 
has been recently applied to study dynamics of weak lo- 
calization in quasi-lD waveguides and 3D slabs 0, and 
dynamics of localized waves in unbounded ID and 2D 
media The key distinctive feature of our approach, 



as opposed to these works, is to allow for position depen- 
dence of the renormalized diffusion coefficient, which is 
of particular importance near open sample boundaries, 
where constructive interferences are strongly suppressed 
due to leaks. Our work may guide experimental quests 
for Anderson localization of ligh t [y, LLJ] and be of help 
in studies of 'random lasers' |2l| . 

In this paper we consider a 3D disordered slab confined 
between the planes z = and z = L 3> L Our theoretical 
model reduces to the diffusion equation for the intensity 
Green's function C: 

[-ifi - V ■ D(r, fi)V] C*(r, r', fi) = 5(r - r') (1) 

with a position- and frequency-dependent diffusivity 



(2) 



where Db = v£/3 and v is the energy transport velocity. 
These equations have to be solved with the boundary 
conditions C^fZo[D (z, {Y)/D b]9 z C = at z — and z = 
L. We adopt zq — 2/3-1 A larger z ~ I allows to account 
for internal reflections at the surfaces of the slab, but does 
not affect our conclusions qualitatively. To regularize the 
unphysical divergence of C(r, r', f2) for r ^ r', we work 
with its Fourier transform C(q_L, z, z' , f2) in the xy plane 
and then introduce an upper momentum cut-off q± = fi/£ 
to obtain C(r, r, fi) in Eq. J5J. The numerical constant 
ji = 1/3 is chosen such that the localization transition in 
the infinite medium occurs at k£ = 1 . 

Following Refs. and 0, we introduce the leakage 
function 

Pt{oi) = — lim T>[z, — ia + e)d z C(0, z, z, — ia + e) 

27T cfO . 



T>(z, —ia — e)<9 z C(0, z, z', —ia — e) 



(3) 



with z = L and z' — £. The Laplace transform of the 
leakage function yields the time-dependent transmission 
coefficient T(t). Deep in the diffuse regime (kl — > oo) 
the leakage function is a sum of delta-peaks located at 
a = n 2 jtn with tr> = (L + 2zq) 2 / it 2 D b the typical time 
needed to cross the disordered sample by diffusion. For 
long times t > to only the first peak is relevant and the 
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FIG. 1: Time-dependent diffusion coefficient for a disordered 
slab of thickness L — 1001 in the regime of weak disorder 
(the product of wavenumber k and mean free path i exceeds 
unity). Dashed lines show the result of the perturbational 
calculation Inset: The first peak of the leakage function 
P T (a). We find P T (a) ~ 6 (a - l/t D ) and D(t)/D B = 1 in 
the limit kl —> oo. 



transmission coefficient Tit) decays exponentially with 
t/t D . 

Let us first see what happens when kl is large but fi- 
nite. We solve Eqs. (0 and (J2J) numerically by extending 
the approach of Ref. [20j to Q ^ and find that the first 
peak of Pt{ol) shifts to smaller values of a while acquir- 
ing a finite width (see the inset of Fig. To characterize 
the deviation of the transmission coefficient from the pure 
exponential decay, it has proven convenient 0,0] to think 
of T(t) as if it were still resulting from a diffusion pro- 
cess though with a time-dependent diffusion coefficient 
D(t) = — [(L + 2z ) 2 /n 2 ](d/d,t)\nT(t). As can be seen 
from Fig. ^ D{t) is smaller than Db and depends on 
time, though its time dependence is rather weak. This is 
also confirmed by the analytic calculation that we per- 
formed assuming t ^> tjj and ki 3> 1 : 
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(4) 



We see that the time-independent deviation of D(t)/DB 
from unity is roughly l/(k£) 2 (in agreement with Rcf. 

where the dependence of V on r was ignored), 
whereas the time-dependent deviation is another fac- 
tor t/L smaller. This implies that in 3D samples non- 
exponential transmission is much more difficult to ob- 
serve experimentally than in quasi- ID waveguides, where 
a linear decay of D(t) with t/tu has been predicted and 
observed in the diffuse regime to < tji 0, 01 ■ The fact 
that D(t) is virtually independent of time in 3D is con- 
sistent with purely exponential T(t) observed in the ex- 



FIG. 2: Time-dependent transmission coefficient T(t) of a 
3D disordered slab of thickness L in the localized regime. For 
each value of ki, curves corresponding to L/£ = 10, 12, and 14 
are shown (£ is the localization length). For clarity, the curves 
corresponding to ki = 0.7 and 0.6 are shifted upwards by 3 
and 6 units, respectively. The t~ (1+a ^ slope (with s « 0.85, 
dashed line) obeyed by all curves crosses over to exponential 
decay for times t > 1/ot. Lower inset: leakage threshold a* 
versus slab thickness L at ki = 0.6 (+), 0.7 (circles) and 0.8 
(triangles). The dashed line is a* = (Db/£ 2 ) exp(— L/£). Up- 
per inset: leakage functions Pr(a) behave as a" (dashed line), 
except for very small a < a* where they vanish rigorously. 



periments of Ref. [£| despite the rather small value of 
ki ~ 4. 

We next turn to the regime of Anderson localization 
(k£ < 1) which is the primary subject of this paper. For 
the parameters used in Fig. |21 we find the peak of T(t) 
at t « to. A detailed study of the arrival time of the 
pulse will be a subject of our future work, but we note 
that this time is determined by relatively short diffusion 
paths, corresponding to large-fi solutions of Eq. (JTJ and 
(J2J), which are only weakly affected by localization ef- 
fects. The decay of T(t) following the peak is power-law: 
T(t) ~ H 1+s ) with s ~ 0.85, in sharp contrast to the 
exponential decay in the diffuse regime. The power-law 
decay of the transmission coefficient originates from the 
power-law growth of the leakage function Pt(o-) ~ ot s 
shown in the upper inset of Fig. [3 The latter is observed 
only for a > a* = {Db/£, 2 ) exp(— L/£) (see the lower 
inset of Fig. 0, where £ = U(kl) 2 /[l - (H) 4 ] < L is 
the localization length. The leakage function Pr(ot) van- 
ishes for a < a*. An accurate analysis yields Pt{o) ~ 
(a — a*) p (with p i=s 0.5) for a j a*, which corresponds 
to T(t) ~ exp(— a*t)/t p+1 for very long times t > I /a*. 
The time-dependent diffusion coefficient 



D(t) 
D B 



(s + \)t D /t, 



t D < t < l/a* 



a*t D + (p+ 1)W*> t>l/a* 



(5) 



decays as 1/t, testifying that the transmission coefficient 
is not pure exponential. The time window tjj <C t < 1 /a* 
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FIG. 3: Leakage function Pr{o) in reflection from a 3D dis- 
ordered half-space in the localized regime. Dashed lines are 
linear fits to the initial parts of the curves. Arrows show ao 
up to which the linear dependence of Pn(a) on a is expected. 
Inset: ao as a function of kl. Dashed lines show approximate 
analytic results a = (L> B /4£ 2 )[1 - 6(kl) 4 ] for kl j and 
a = (8D S /27£ 2 )(1 - kl) for kl T 1. 



FIG. 4: Distribution of leakage rates P(T) for a slab of thick- 
ness L in the localized regime (localization length f). For 
clarity, only the results for kl = 0.8 are presented. For other 
kl < 1 the curves are identical, except that the region of 
diffuse behavior Db/3£ 2 < T < Db/1 2 shrinks when kl de- 
creases and £ approaches I. 



of power-law decay of T(t) opens up progressively as £ 
decreases below L. The localization transition at kl = 1 
is smooth, in contrast to the infinite medium where kl = 
1 defines a sharp boundary between diffuse and localized 
regimes of wave propagation. 

Let us now consider localization effects in time- 
dependent reflection, first focusing on a disordered half- 
space {L — > oo). For diffuse waves {kl > 1), the aver- 
age reflection coefficient R(t) of the half-space scales as 
1/t 3 / 2 . Such a scaling of R(t) corresponds to the leak- 
age function Pn(a) ~ y/a, where the leakage function 
in reflection Pr(o) is defined by Eq. © with z = 
and with a change of sign. The square root scaling of 
Pr{oi) is also present in the localized regime, but only 
for a > a ~ D B /£, 2 , yielding R(t) ~ f/< 3/2 for short 
times t <C 1/ao. The waves leaving the sample at such 
short times have not penetrated deeply enough to be af- 
fected by localization effects and it is therefore natural 
to recover the diffusion result. In contrast, Pr(o) scales 
linearly with a for a < a (see Fig. |3J), leading to the 
scaling law R(t) ~ 1/t 2 for long times t ^> 1/ao- Inter- 
estingly, the value of ao can be calculated analytically, 
although the result is quite lengthy (see the inset of Fig. 

EJ). 

The inverse quadratic scaling of the reflection coeffi- 
cient with time was reported first for ID semi-infinite 
disordered media j^- The transition from R(t) ~ 1/t 3 / 2 
to R(t) ~ 1/t 2 also occurs in quasi-lD disordered waveg- 
uides |j,|23j. The 1/t 2 scaling of the reflection coefficient 
seems therefore to be a hallmark of Anderson localiza- 
tion. For a medium of finite thickness L>( this scaling 
of R(t) should still be valid until a large time 1/a*, be- 



yond which an exponential decay sets in. 

An alternative way of dealing with Anderson localiza- 
tion is to study (quasi-)modes of the disordered sam- 
ple. The statistical distribution P(T) of decay rates T 
of the modes has recently attracted considerable atten- 
tion I24I Eil EH. I27I l2sj , also in the context of random 
lasing 29]. In the localized regime, the Laplace trans- 
form of P(T) approximately equals the so-called sur- 
vival probability P s (t) The latter, in turn, obeys 
dP s (t)/dt — —2T(t), where T(t) is the Laplace transform 
of the leakage function Pt{&) given by Eq. averaged 
over the source position z' inside the disordered sample. 
This yields P(T) = 2P T (T)/T. The distribution of leak- 
age rates estimated in this way is (see Fig. 

P(T) = (6) 

0, r < a* 

K/L) x 1/r, a* < r < D B /se 

(i/3)V^xi/r 3 / 2 , D B /3e < r < D B /l 2 

The 1/r behavior of P(T) is a hallmark of Anderson 
localization and has already been report ed p reviously for 
both disordered and chaotic systems [24. l28j (note that a 
slightly different result P(T) ~ 1/r 1,25 has been reported 
by Titov and Fyodorov 25]). P(T) ~ _l/r 3 / 2 is typical 
for diffuse regime of wave propagation |2fJ, |28j . This part 
of P(r) is due to quasi-modes located closer than £ to 
the boundaries of the sample, and which therefore behave 
more as extended than as localized modes. T > D B /l 2 
correspond to modes with lifetimes shorter than the mean 
free time and cannot be correctly described by our theo- 
retical model. 

The mode picture of Anderson localization allows a 
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better understanding of 1/t 2 scaling of short-pulse re- 
flection from a disordered half-space. A short pulse at 
z = 0, t = excites exponentially localized modes of dis- 
ordered sample with relative weights exp(— z/£), where z 
is the position of a given localized mode. Each mode then 
leaks exponentially with rate T proportional to its inten- 
sity at the surface of the medium 

Hence, the relative weight of a given mode is proportional 
to its leakage rate T, which immediately yields Pr(o) ~ a 
and R(t) ~ 1/t 2 , if interferences of different modes are 
ignored. 

According to our theory Pr,t{oi) and P(T) all ex- 
hibit gaps for < a, T < a*. This seems to contra- 
dict the arguments of Refs. H 0, 0, HI] that predict 
P(T) - exp[- ln d (Tt H )] (with d = 2 or 3) in the limit of 
small r, due to some special, rare realizations of disorder. 
However, we consider a medium with infinite Heisenberg 
time and therefore our results cannot be directly com- 
pared to those of Refs. HUES El]. The unification of 
smalhr results corresponding to tjj — ► oo and tn < oo 
constitutes a major challenge for future research. 

In conclusion, we have presented a self-consistent the- 
oretical approach to the dynamics of Anderson localiza- 
tion in open 3D disordered media. We have calculated 
the time-dependent transmission and reflection, as well 
as the distribution of mode decay rates of open disordered 
slabs in the localized regime. The transmission exhibits 
power-law decay in time, followed by an exponential de- 
cay at very large times. The reflection crosses over from 
1/t 3 / 2 to 1/t 2 behavior. The distribution of mode de- 
cay rates P(T) has essentially two parts: P(T) ~ 1/r 
and P(T) ~ 1/r 3 / 2 , corresponding to modes localized 
deep inside the sample and near the boundaries, respec- 
tively. These results are of particular interest in the con- 
text of recent experiments on Anderson localization of 

light ism. 
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